Chiral d-wave RVB state on honeycomb lattice as a generalized staggered flux phase 



O 
(N 



I 

o 
o 



> 
(N 



Tao Li 

Department of Physics, Renmin University of China, Beijing 100872, P.R.China 

(Dated: July 20, 2011) 

We show the chiral d-wave RVB state on honeycomb lattice stands as a natural generalization 
of the staggered flux phase on square lattice. Although the state is generated from a time reversal 
symmetry broken mean fleld ansatz, it actually represents a fully symmetric spin liquid state with a 
positive definite wave function in the sense of Marshall sign rule for unfrustrated antiferromagnets. 
The evolution of the state with the parameter A/x follows exactly the same manner as that of the 
staggered flux phase on square lattice. The critical pairing strength corresponding to the 7r-flux 
phase is found to be A/x = \/2. As a result of the geometric frustration between neighboring 
plaquette on honeycomb lattice, a direct generalization of the U{1) staggered flux pattern on square 
lattice to honeycomb lattice is impossible. Replacing it is the chiral d-wave state with Z2 gauge 
structure. However, this Z2 gauge structure is found to be ineffective after Gutzwiller projection 
and the system does not support topological degeneracy. The chiral d-wave RVB state is also found 
to be a rather good variational state for the Heisenberg model on honeycomb lattice. The spin 
correlation of the chiral d-wave state is found to be greatly enhanced as compared to the mean fleld 
prediction. 

PACS numbers: 75.10.Kt,71.27.-|-a 



I. INTRODUCTION 

The study of correlated electrons on honeycomb lattice 
has attracted much interest recently. The honeycomb lat- 
tice has the smallest possible coordinate number of 3 for 
a two dimensional lattice and hosts Dirac-type disper- 
sion at half filling. It is interesting to understand how 
the electron correlation effect will manifest itself in such 
a background. 

A comparison with the square lattice case is of par- 
ticular interest. Both lattices are bipartite and the an- 
tiferromagnetic long range order is unfrustrated in the 
strongly correlated limit. However, at half filling, the 
honeycomb lattice has a Dirac-type dispersion and a fi- 
nite correlation strength is needed to induce antifcrro- 
magnetic order on honeycomb lattice, while on square 
lattice, the system form antiferromagnetic long range or- 
der immediately as we turn on interaction as the half 
filled system has a nested Fermi surface. More recently, 
it is discovered that the transition from the semimetal 
phase to the antiferromagnetic ordered state on honey- 
comb lattice may take place in a two step manner, leaving 
a small but finite intermediate correlation range in which 
the system may enter a spin liquid state with unconven- 
tional characteristics 

On square lattice, the antiferromagnetic order is de- 
stroyed by charge carrier doping. The resultant state 
is believed to have d-wave superconducting pairing and 
may be responsible for the superconductivity of the high- 
Tc cuprates. What is the analogy of these doping effect 
for the correlated electron system on honeycomb lattice? 
More specifically, what is the pairing symmetry of the re- 
sultant superconducting state, provided that the doped 
system does enter some kind of unconventional super- 
conducting state. This inquire is even more attractive on 



honeycomb lattice than on square lattice, as the undoped 
honeycomb system already show signature of spin liquid 
ground state, which is nothing but a undoped version of 
the pairing state, the RVB state. 

Superconductivity on system with a honeycomb lat- 
tice, especially that on graphene, has been studied by 
many researchers from both theoretical and experimen- 
tal perspectives [il-Q. However, in most of these studies, 
the electron correlation effect does not play an essential 
role and the related pairing state has either s-wave or 
extended s-wave pairing symmetry, originated from ex- 
changing phonon or other excitations. When the electron 
correlation effect is taken into account, the chiral d-wave 
pairing stated, been proposed as a generalization 
of the d-wave pairing state on square lattice to honey- 
comb lattice(see Fig.l). However, it is not clear how the 
time reversal symmetry broken chiral state are connected 
to the antiferromagnetic ordered state at half filling. 

At this point, it is quite useful to recall how the pair- 
ing symmetry of the superconducting state is determined 
in the strong correlation limit on square lattice. Here, a 
RVB state called staggered flux phase is found to be the 
best variational state at half filling in the restricted space 
of Fermionic RVB state 10- 13] , although it is known that 
the half filled system inevitably has an antiferromag- 
netic long range order. The staggered fiux phase has 
the virtue that the antiferromagnetic short range corre- 
lation is greatly improved by the introduction of a non 
zero staggered flux. Up to a SU(2) gauge transforma- 
tion, the introduction of the staggered flux is found to be 
equivalent to the introduction of d-wave pairing between 
the Fermionic spinous, especially, when the pairing and 
hoping order parameter becomes identical, the staggered 
flux becomes tt. For general value of the staggered flux, 
the pairing amplitude and the flux is related by the rela- 
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tion tan 4 = • It is thus quite natural that the doped 
system will choose the d-wave pairing state as its super- 
conducting state as the antiferromagnetic superexchange 
is already optimized by the staggered flux. 

Thus, a natural questions arises as what is the ana- 
log of the staggered flux phase on the honeycomb lattice, 
which can hopefully improve the superexchange interac- 
tion by some kind of pairing(flux). However, except for 
the TT-flux phase, a naive generalization the staggered flux 
pattern for general flux value to the honeycomb lattice 
is impossible as a result of the geometric frustration be- 
tween the neighboring hexagonal plaquette(see Fig. 2). 
Such a difference with square lattice originates from the 
frustrated nature of the dual lattice of the honeycomb 
lattice, which is the triangular lattice. Thus a general- 
ization of the staggered flux phase to honeycomb lattice 
in the space of U{1) staggered flux is impossible. How- 
ever, generalization of the staggered flux phase in the 
larger space of SU (2) flux is still possible. In this paper, 
we show that the RVB state with chiral d-wave pairing 
just fit such a need. 

More specifically, we show that the chiral d-wave RVB 
state on honeycomb lattice can be taken as a general- 
ized staggered flux phase. It actually describes a spin 
liquid state with the full symmetry of the Heisenberg 
model on honeycomb lattice and evolves in exactly the 
same way as the staggered flux phase on square lattice 
as a function of the parameter A/x(see Fig. 3). When 
— = 0, the state reduces to the uniform RVB state on 
honeycomb lattice, which has Dirac-type spinon disper- 
sion. When - = \/2, the state is gauge equivalent to the 
TT flux phase on the honeycomb lattice. When exceed 

•\/2 and approaches infinite, the state again evolves back 
to the uniform RVB state with zero fiux. It is interesting 
to note that in an earlier work on symmetric spin liquid 
state on honeycomb lattice, the chiral d-wave RVB state 
is classified to be a Z2 spin liquid in the neighborhood of 
the uniform RVB state [14j. As its cousin on square lat- 
tice, the introduction of the staggered fiux improves the 
antiferromagnetic correlation and we find the chiral d- 
wave RVB state stands as a rather good variational state 
for the Heisenberg model on honeycomb lattice. Unlike 
the staggered flux phase on square lattice, for general flux 
value, the mean field ansatz for the chiral d-wave RVB 
state posses a Z2 rather than U{1) gauge structure. We 
also find that as a result of the bipartite nature of the 
lattice, both the staggered fiux phase on square lattice 
and the chiral d-wave RVB state on honeycomb lattice 
satisfy the Marshall sign rule for unfrustrated antiferro- 
magnet. Such sign structure will prohibit the system to 
show topological degeneracy in the Z2 sense, no matter 
what is the gauge structure of the mean field ansatz for 
the RVB state. 

This paper is organized as follows. In the next two 
sections, we will review some well known results about 



the RVB state in general and the staggered flux phase 
on square lattice in particular to set up the stage for our 
discussion. In section IV, we present our results for the 
chiral d-wave RVB state. Section V contains a detailed 
comparison between the staggered flux phase on square 
lattice and the chiral d-wave RVB state on honeycomb 
and a discussion of some related issues. 




FIG. 1: The phases of the pairing potentials in the chiral 
d-wave RVB state on honeycomb lattice. 




FIG. 2: The geometric frustration of (7(1) staggered flux on 
the honeycomb lattice. 



II. RVB STATES AND THEIR GAUGE 
STRUCTURES 

The RVB states studied in this paper are those de- 
rived form Gutzwiller projection of BCS-type mean 
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FIG. 3: The evolution of the gauge structure for the staggered 
flux phase on square lattice(a) and the chiral d-wave RVB 
state on honeycomb lattice(b). Here GW-Dirac denotes the 
uniform RVB state on honeycomb lattice, which has a Dirac- 
type spinon dispersion. 



field ground state of the following general mean field 
ansatz[l33 
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(1) 
Eq.(l) 



can also be taken as the saddle point approximation of 
a field theoretical description of the RVB state, in which 
Xi,j and are interpreted as two types of RVB or- 
der parameters. In the field theoretical formulation, the 
Gutzwiller projection into the subspace of no double oc- 
cupancy amounts to integrating over the time component 
of the gauge fluctuation. 

It is important to note that within the subspace of 
no double occupancy, a description of the RVB state in 
terms of Xi,j ^-nd A^j- becomes redundant as a result of 
the gauge degree of freedom in the mean field ansatz . 
Such redundancy is accompanied by an underlying gauge 
structure of the RVB state so constructed. To see more 
clearly, it is better to rewrite the mean field ansatz in the 
Nambu form 



H 



ihj) 



^3 



(2) 



in which -0^ 



A* 



(ci_-^,c|^) is a two component spinor. 
is a 2 X 2 matrix. It is clear then 



i,j Xi,j 

that the system is invariant under the following SU{2) 
gauge transformation 



(3) 



2 J ' 



in which Wi is a site-dependent SU (2) matrix. As only 
the SU{2) gauge singlet survive the Gutzwiller projec- 
tion, Uij and WiUijWj actually describes the same RVB 
state. 

Although Uij is not gauge invariant and thus unphys- 
ical, it does contains important gauge invariant informa- 
tion on the structure of the RVB state. To uncover such 
internal gauge structure of the RVB state, it is a common 
practice to construct the loop operators [l^, Pi{C) = 
UijUj^kUk,i ■ ■ ■ Um,i, in which i denotes the starting point 
of the loop C and j,k, - ■ - m are the remaining sites along 
the loop. Under the SU{2) gauge transformation, a loop 
operator transform as Pi{C) — WiPi{C)W- and thus its 
trace form a gauge invariant quantity. A RVB state is 
called to have SU{2) gauge structure if all loop operator 

are proportional to the identity matrix ''o = ^ q ]' 

Otherwise, it is called to have a [/(I) gauge structure if 
all Pi{C) starting from any given site i commute with 
each other. If there are more than two loop operators 
starting from the same site do not commute with each 
other, the RVB state is called to have Z2 gauge struc- 
ture. The gauge structure of the RVB state is closely 
related to the existence of soft gauge mode in the long 
wave length limit in the effective field theoretic descrip- 
tion of the RVB state. For example, if the gauge struc- 
ture of the RVB state is SU{2), then there will be SU{2) 
soft gauge mode in the long wave length limit. While if 
the gauge structure of RVB state is Z2, then there is no 
soft gauge mode at all in the long wave length limit. In 
such a case, the predictions made at the level of saddle 
point approximation is believed to be more reliable than 
those cases in which the RVB state has a U{1) or SU{2) 
gauge structure. In addition, the system with Z2 gauge 
structure is argued in the effective field theory to support 
topological order and topological degeneracy - existence 
of degenerate ground states that can not be differentiated 
from any local probe [l5|. 

The RVB state is generated from Gutzwiller projection 
of the BCS-type mean field ground state 



|RVB) = PgIBCS). 



(4) 



Such a state can in general be rewritten in the form of 
condensed spin singlet pairs 



|RVB)=PG^a(*-j)44J |0), 



(5) 
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in which a{i — j) is the so called RVB amplitude and 
satisfy the relation a{i — j) = a{j — i), N is the number 
of lattice site. As gauge equivalent mean field ansatz Uij 
lead to identical RVB state, a RVB state is said to be 
symmetric if and only if the new ansatz after symmetry 
transformation is gauge equivalent to the old one. The 
projective symmetry group(PSG) provides a convenient 
scheme to classify the symmetric RVB state 15 [ . 

For RVB state with Z2 gauge structure, which is ar- 
gued to posses topological degeneracy from effective field 
theory, a Z2 topological excitation called vison can be 
constructed by reversing the sign of RVB amplitude on 
bonds that crossing a branch cut line originated from the 
center of the vison an odd number of times. The topo- 
logical nature of the vison excitation can be seen from 
the fact it acts as a tt flux tube located at its origin for 
spinon traveling around it. Especially, when the RVB 
state is defined on a torus and a vison is trapped in one 
of the holes of the torus, one is left with a state that 
globally distinct from, but locally indistinguishable(and 
thus degenerate in the thermodynamic limit) from the 
original RVB state. This is nothing but the topological 
degeneracy. However, on bipartite lattice, the topologi- 
cal degeneracy predicted by the effective field theory can 
be lifted as a result of some special phase structure on 
the RVB state. 

The detection of the Z2 topological degeneracy is sim- 
plified by the fact that trapping a vison in the holes of the 
torus is gauge equivalent to change to boundary condition 
from periodic to anti-periodic or vice versa in the mean 
field ansatz 16|, ll7| . Thus to see if a particular RVB state 



does posses topological degeneracy in practice, we only 
need to calculate the overlap between RVB states gen- 
erated from mean field ansatz with different boundary 
conditions around the holes of the torus and extrapolate 
the result to the thermodynamic limit. If the overlap ex- 
trapolate to zero in the thermodynamic limit, then the 
RVB state is said to posses topological degeneracy. 



III. THE STAGGERED FLUX PHASE ON 
SQUARE LATTICE 

The staggered flux phase (or d-wave RVB state) on 
square lattice plays a very important role in our under- 
standing of the high-Tc superconductivity in cuprates. 
In the following we will review briefly some of the most 
important properties of this state [l3| for a comparison 
with the chiral d-wave RVB state on honeycomb lattice. 

The staggered flux phase on square lattice is generated 
from the following mean field ansatz 

HsF- - E (e^^'-4,c„.-t-/i.c.), (6) 

in which the phase factor is introduced to guarantee 
that each plaquette of the square lattice is threaded by 



a U{1) flux of value ±$ arranged in a staggered pattern. 
Although the mean field ansatz breaks the time rever- 
sal symmetry for general value of it is well known 
that the staggered flux phase is gauge equivalent to the 
d-wave RVB state generated from the following time re- 
versal symmetric d-wave BCS mean field ansatz 



Hd-Bcs = -X ^4,a(^j,'y + h-c.) 



(7) 



in which A^j- = ±A for nearest neighboring sites i and j 
in the x/y direction. Here both x and A are real and is 
determined from $ by 

tan- = -. 8 

4 X 

In the following, we will use the d-wave gauge. 



A. Mean field description 

The spinon excitation spectrum in the staggered fiux 
phase has the Dirac-type linear dispersion and is given 
by Ek = v/(ek)^ + (Ak)^, in which ek = -2x(cosk::c -I- 
coskj^), Ak 2A(cosk2; — coskj,). The nodes are located 
at (ka;,kj^) = (±7r/2, ±7r/2). As we will be concerned 
with the spin structure factor of the RVB state below, 
we present the mean field prediction for it here. The 
spin structure factor is defined as 



(9) 



The spin structure factor in a BCS mean field ground 
state is given by 



Siq) 



8iV2 ^ 

k 



(1 



CkEk+q + AkAk^ 
E-kEk+a 



(10) 



Since the square lattice is bipartite and ek = — ek+q, 
Ak = -Ak+Q for Q = (7r,7r), we find S'(Q) = ^ in the 
mean field theory. The spin structure factor is thus in- 
dependent of the value of the staggered flux in the mean 
field theory. It also extrapolates to zero in the thermo- 
dynamic limit. 



B. Gauge structure 

The loop operator around the elementary plaquette of 
the square lattice in the staggered flux phase is given by 

Pl{C) = U^UyU^Uy 

= {x'-A'r^4x'A'+4^xA{x'-A')r2, 
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FIG. 4: The loop operator defined on the elementary plaque- 
tte of square lattice. 



in which Ux = -X^a + ^'^i. Uy = -x^i, - An. Here ti, 
T2 and ra denote the three Pauh matrixes in the internal 
SU{2) space. 

When A/x = 0, |A| = |x| or x/A = 0, the loop op- 
erator is proportional to To(here we assume that x and 
A satisiy the normalization condition x^ + A^ = 1). By 
induction we can show that all loop operators are propor- 
tional to To at these three points. The state with x/A = 
and that with A/x = are in fact gauge equivalent and 
is the well known uniform RVB state on square lattice. 
The state with |A| — |x| is the so called 7r-flux phase on 
square lattice. Thus, both the uniform RVB state and the 
TT-flux phase posses a SU{2) gauge structure. For gen- 
eral value of the staggered flux, the loop operator is not 
proportional to tq and gauge structure can be shown to 
be C/(l). The evolution of the gauge structure with A/x 
in the staggered flux phase is summarized in Fig. 3(a). 



staggered flux $ and — As the gauge flux is defined 
modula 27r, the state with staggered flux $ and 27r— $ are 
gauge equivalent. In the d-wave gauge, this indicates that 
the state with A/x = a and the state with A/x = 1/a 
are gauge equivalent. 



D. Sign structure of the staggered flux phase and 
topological degeneracy 

In an Ising basis for the spins, the staggered flux phase 
can be expanded as follows 



|SF) = ^V^({aJ)|{aJ), 



(11) 



in which |{cri}) — |cri, • • • , (Tat) denotes an Ising basis. 
Then the wave function ^{{cTi}) can be shown to be 
real and satisfy the Marshall sign rule 17 , 1^ for unfrus- 
trated antiferromagnet, namely, the sign of is given by 
(—1)^ ^, with N^^ denoting the number of down spins 
in A sublattice. 

Although the mean field ansatz of the staggered flux 
phase has at least a U{1) gauge structure, it is still pos- 
sible to construct a state with a trapped vison in the 
holes of a torus and check its orthogonality with the state 
with no trapped vison. Such a calculation has been done 
in previous studies fio'l and it was found that the over- 
lap extrapolates to a finite value in the thermodynamic 
limit, indicating no topological degeneracy for the stag- 
gered flux phase. This result is claimed to be a strong 
support of the effective field theory argument. However, 
a later investigation [l^ indicates that the Marshall sign 
rule satisfied by the staggered flux phase plays a more 
essential role in removing the topological degeneracy. 



E. Variational energy and spin structure factor 



C. Symmetries 

The staggered flux actually respects all the physical 
symmetries of the Heisenberg model on square lattice. 
While translational symmetry and inversion symmetry is 
manifest in the d-wave gauge, the four-fold lattice rota- 
tional symmetry can be shown as follows. 

Under the four-fold lattice rotation, the mean field 
ansatz transform as follows: x ~^ X: A — —A. A 
global SU{2) gauge transformation with Wj = ir^ suf- 
fices to recover the original ansatz. The rotated ansatz 
and the original ansatz thus describe the same RVB state. 
This prove the rotational symmetry of the staggered flux 
phase. 

The four-fold rotational symmetry of the staggered flux 
also implies the gauge equivalence between the state with 



Both the short range and the long range spin correla- 
tion are seen to reach extreme at the uniform RVB state 
and the 7r-flux phase. This is a result of the gauge equiv- 
alence between states with staggered flux $ and — $ and 
that between states with staggered flux $ and 27r — 
from which one can easily show that the spin correlation 
function should reach their extreme in the uniform RVB 
state and the tt-Aux phase. 

The correlation between spins on nearest neighboring 
sites, S(i j) i'^i ' "^i) 1 which is proportional to the vari- 
ational energy for the Heisenberg model on square lattice, 
is presented in Fig. 5 as a function of A/x. The optimal 
value for A/x is found to be about 0.3. The variational 
energy reaches local maximum in both the uniform RVB 
state and the tt-Aux phase. 

To probe the long range spin correlation, we have cal- 
culated the spin structure factor of the staggered flux 
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FIG. 5: The variational energy per bond calculated for the 
Heisenberg model on square lattice as a function of A/x- 
The calculation is done on a 16 x 16 lattice with periodic- 
antiperiodic boundary condition. The error bars are smaller 
than the size of the symbols. 
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FIG. 6: The spin structure factor at q = (tt, tt) as a function 
of A/x for the staggered flux phase on square lattice. The cal- 
culation is done on a 16 x 16 lattice with periodic-antiperiodic 
boundary condition. The error bars are smaller than the size 
of the symbols. The red line denotes the mean field prediction, 
which is independent of A/x- 



phase. The spin structure factor peaks at the antifer- 
romagnetic ordering wave vector Q — (tt, tt). In Fig. 6, 
we present the result for S{Q) as a function of A/x. 
The spin structure factor decrease monotonically with 
the staggered flux for $ < tt. The small peak close to 
the uniform RVB state is a finite size effect and it moves 
toward the uniform RVB state with increasing lattice size 
and disappears in the thermodynamic limit. 

In a recent work [2^, we have shown that the uniform 
RVB state on square lattice actually describe a state 
with antiferromagnetic long range order as a result of the 
nested spinon Fermi surface. The spin structure factor 
follows the S{Q) w S'o + a(l/i)^/* behavior as a function 
of the linear scale of the lattice, L. While in the 7r-flux 
phase, S{Q) is found to follow the 5(Q) « be- 
havior and no magnetic order is detected. In between 
the two extremes, the staggered flux has no magnetic 
long range order and the spin correlation function decay 
algebraically with distance with an exponent depending 
on the value of the staggered flux. Thus the uniform RVB 
state state on square lattice is an isolated singular point. 



IV. THE CHIRAL D-WAVE RVB STATE ON 
HONEYCOMB LATTICE 




FIG. 7: The unit cell of the honeycomb lattice. 



in which denotes sum over nearest neighboring 

sites on honeycomb lattice. Aij denotes the pairing order 
parameter between site i and j. In the chiral d-wave RVB 
state, Aij is a complex number and the phase of Aij in 
the three directions differ with each other by ±^(see 
Fig.l). 



The chiral d-wave RVB state on honeycomb lattice is 
generated from the following BCS mean held ansatz 



A. Mean field description 



Hbcs ^-t i4,aCj,.+h.c.)+ J2 ^^ji'^clAs+h-C-), 

The mean field Hamiltonian for the chiral d-wave state 
(12) can be easily diagonalized in momentum space in which 
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it takes the form 



Hbcs ^Jl^k^^i' 



(13) 



in which -0^ is a four component spinor given by V'jf 



The 4x4 matrix M is given by 



Here A and B are sublattice index. 



M : 



/ 7k Ak \ 

A_k 

Alk -7k 

-7^ / 



7k 




(14) 



in which 7k = (1 + e''^-^^ + e"'''-'^^), Ak = A(l + 
gi--r g-^'^-Si -(- e*"5^e"*'^ '^^). Here (5i and (52 are the two 
lattice translational vectors of the honeycomb lattice(see 
Fig.7). 

The eigenvalues of the Hamiltonian are given by 




FIG. 8: The loop operators of the chiral d-wave RVB state on 
the elementary plaquette of the honeycomb lattice. The gray 
dots denotes the common starting point of the loop operators. 



Ek± 



7^ 



|A_ 



±^/W^ 



12^,2 



)2+4|5k|^ 



(15) 



7-k^ 



7kAj: 



in which gk 

When A = 0, the dispersion reduces to that of 
the Dirac Fermion, Ek = ±|7k|- At (^1,^2) = 
±(27r/3, — 27r/3), Ek — 0. Here we have used the con- 
vention that k — fcibi + fc2b2 for the momentum, bi 
and b2 are the two reciprocal vectors. They are dual to 
the lattice translational vectors 61 and 62 and satisfy the 
relation Si ■ hj — 6ij. 

When (/ci,fc2) = ±(27r/3, -27r/3), 7k = and thus 
g{k) — 0. It can also be shown that one of the 
two gap functions, Ak or A_k, is zero at {ki,k2) — 
±(27r/3, — 27r/3). It is then easy to verify that the nodes 
at {ki,k2) — ±(27r/3, — 27r/3) persist for any value of A. 
In addition, new nodes may emerge at nonzero A. In the 
appendix, we show the dispersion for some typical values 
of A/x. It should be noted that the location of the addi- 
tional nodes move continuously in the momentum space 
with A/x- 

On the honeycomb lattice, the spin structure factor at 
the antiferromagnetic ordering wave vector can be de- 
fined as follows 



(16) 



B. Gauge structure 

From every sites of the honeycomb lattice, there are 
three elementary loop operators. They are given by 



P{Ci) = {UiU2U3f 



(17) 



PiC2) 



{U2U3Uif 
{U3UiU2f, 



in which Ui 



VIA _ 
2 '2 



'XT3 + Ari = pui, U2 = -XT3 - yTi + 

pu2, U3 = ~XT3 - - ^^T2 = pu3 are 
RVB order parameters on bonds in three directions, p — 
\/x^ + A^ is the 'length' of the RVB order parameter. 
Mi,M2,U3 are three Hermitian matrix satisfying uiUi — 
U2U2 = U3U3 — tq. The products of Ui, U2 and U3 in 
Eq.l7 are given by 



U1U2U3 



U2U3U1 = 



.3^3 .2 



(18) 



A2 



){-XT3 + Ati - V3AT2) 



.3\/3 



XA' 



ix' 



-){-XT3 + An + V3AT2) 



U3U1U2 = -I 



.3V3 



XA^ 



From the mean field theory, it is straightforward (though 
somewhat tedious) to show that S{Q) is independent of 
A/x and is given by the S{Q) = As for the staggered 
flux phase on square lattice, the A/x independence of 
S{Q) on honeycomb lattice at the mean field level is a 
result of the bipartite nature of the lattice. 



ix' 



-){-XT3-2An). 



From Eq.l8, it can be easily seen that at the three 
special points A = 0, A = "\/2x and x = 0(with the nor- 



malization condition A^ - 



■X 



1 assumed) , all the three 



loop operators are proportional to the identity matrix tq. 



8 



Sublattice/Parity 


(0,0) 


(o,e) 


(e,o) 


(e,e) 


A 


To 


— ms 


—iU2 


—mi 


B 


UlT-s 


— U-2T3 




in 



TABLE I: The gauge transformation that relate the chiral d- 
wave RVB state with A/x = \/2 and the 7r-fiux phase shown 
in Fig.9 



By induction it can then be shown that all loop opera- 
tors are proportional to tq. The system thus has SU{2) 
gauge structure at these three special points. This is ex- 
actly what happens in the staggered flux phase on square 
lattice. In fact, it can further shown that the chiral RVB 
state on honeycomb lattice with A = (x = 0) and 
A = V2x are just the uniform RVB state and tt-Aux 
phase on honeycomb lattice. 

The gauge equivalence between the state with A = 
and X = 1 {Ui = U2 = U3 = —T3) and that with 



and A 



1 (t/i 



X 

U3 — — 5T1 — -3^12) can be shown by the following SU{2) 
gauge transformation. First, the phase of the pairing or- 
der parameter in the state with x — and A = 1 can be 
gauged away by a U{1) gauge transformation. After this 
U{1) gauge transformation, the mean field ansatz takes 
the form Ui — U2 = U3 = ti. Then a uniform SU{2) ro- 
tation with Wi = e^*~'^^ suffices to transform the ansatz 
into the form Ui = U2 = U3 — —T3. Similarly, the gauge 
equivalence between the state with A — V^x and the 
TT-fiux phase on honeycomb lattice(see Fig.9) can be es- 
tablished with the SU (2) gauge transformation given in 
Table. 1 . 




tion in the dual lattice of honeycomb lattice. To see the 
evolution of the gauge structure more clearly, we rewrite 
the loop operator in the following form 

(19) 



P(Ci) = 




PiC2) - 









in which the absolute value of the gauge fiux $ is deter- 
mined by 



cos ■ 



4> _ 3V3xA2 
2p3 ■ 



(20) 



fti, n2 and ^3 are three vectors of unit length and are 
given by 



ni = r{A,-V3A,~x) 
n2 = r(A,V3A,-x) 



(21) 



r(-2A,0,-x), 



in which r = 1/ \/x'^ + 4A2. The three vectors form the 
equal angle with each other as shown schematically in 
Fig. 10. When A = 0, all the three vectors lies in the 
— Ta direction. With increasing A/x, the three vectors 
split from each other. When A/x —> 00, ni, n2 and 713 
evolves into the ti — T2 plane with a 120° angle between 
each other. 















^1 













FIG. 9: The mean field ansatz for the 7r-flux on honeycomb 
lattice. The bonds in red have a negative hopping integral, 
while other bonds have a positive hopping integral. 

For other value of A/x, the three loop operators do 
not commute with each other and the system has a Z2 
gauge structure. This is quite different from the square 
lattice where the SU (2) points are connected by interme- 
diate state with U{1) gauge structure. The origin of such 
a difference can be traced back to the geometric frustra- 



FIG. 10: The direction of the three loop operators in the 
internal 5(7(2) gauge space. 

The structure of the loop operators presented above 
can be viewed as a generalized staggered fiux pattern on 
honeycomb lattice. The absolute value of the fiux for the 
staggered flux phase on square lattice and the chiral d- 
wave RVB state on honeycomb lattice are both uniform. 
It is the direction of the flux that is responsible for its 
staggered character. On square lattice, where the dual 
lattice is also square lattice, a two sublattice (antiferro- 
magnetic) arrangement of the direction of the fiux is the 
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most natural choice. While on honeycomb lattice, where 
the dual lattice is triangular lattice, a three sublattice 
arrangement of the direction for the flux is more natural. 

The discussion above on the gauge structure of the 
chiral d-wave RVB state is summarized in Fig. 3b. Its 
similarity with the figure for the staggered flux phase on 
square lattice is apparent. It is for this reason that we 
call the chiral d-wave RVB state a generalized staggered 
flux phase on honeycomb lattice. 



C. Symmetries 

As will be clear below, the term 'chiral' is in fact not 
quite accurate. It can be shown that the time reversal 
symmetry breaking manifested in the mean field ansazt is 
in fact an artifact of the mean field description. In fact, it 
can be more generally proved that the chiral d-wave RVB 
state actually respects all the physical symmetries for the 
Heisenberg model on honeycomb lattice. The effect of the 
symmetry transformation on the mean field ansatz is to 
induce permutations between Ui, U2 and U3. Thus to 
prove the symmetry of the chiral d-wave RVB state it 
suffice to show that the RVB order parameters after such 
permutation are gauge equivalent to the original one. 

For example, under the time reversal transformation, 
U2 and C/3 are exchanged. Such a change can also be 
induced by the following two-step gauge transformation. 
First, we rotate the ansatz uniformly along the T2 axis 
by TT. Then as the lattice is bipartite, we effect a C/(l) 
gauge transformation by tt uniformly on sites in the B 
sublattice. This proves the time reversal symmetry of 
the chiral d-wave RVB state. Note that the bipartite 
nature of the lattice is essential for the restoration of the 
time reversal symmetry. Another example is provided by 
the 2tt/3 rotation of the lattice, which induces a cyclic 
exchange among Ui, U2 and C/3. It can be easily shown 
that such a cyclic exchange can also be induced by a 
uniform rotation of the ansatz along the Ta-axis in the 
SU{2) space by an angle 27r/3. This proves the rotational 
symmetry of the chiral d-wave RVB state. 

We note that as the exchange between U2 and U3 and 
the cyclic exchange among Ui,U2 and U3 exhausts the 
generators of the permutation group of Ui,U2 and C/3, 
we have in fact proved the full symmetry of the chiral 
d-wave RVB state. 



D. Sign structure of the chiral d-wave RVB state 
and topological degeneracy 

Beside being time reversal symmetric, it can also be 
shown that the chiral d-wave RVB state on honeycomb 
lattice actually has a positive definite wave function in 
the sense of Marshall sign rule for bipartite antiferromag- 
net, although the mean field ground state break the time 




FIG. 11: The overlap between |^'io) and I'I'oi) on a L x L x 2 
honeycomb lattice as a function of the lattice size L. The 
error bars are smaller than the size of the symbols. The oscil- 
lation of the curve is caused by the complex nodal structure 
of the chiral d-wave RVB state, which move continuously in 
the momentum space with A/x- 



reversal symmetry and has a complex valued wave func- 
tion. This remarkable result is in fact a general property 
for all RVB states generated from bipartite mean field 
ansatz 



17f ■ It is also argued that the Marshall sign struc- 



ture will remove from the projected wave function the 
topological degeneracy, even if the mean field ansatz has 
a Z2 gauge structure. 

The absence of the topological degeneracy on bipartite 
system can be argued as follows. As changing the bound- 
ary condition around the hole of a tours will not change 
the bipartite nature of the mean field ansatz, both the 
projected state with trapped vison and without trapped 
vison will be positive definite in the sense of Marshall sign 
rule. It is just this insensitivity of the sign structure to 
the trapped Z2 gauge flux that is responsible for the ab- 
sence of the topological degeneracy. In a recent work[2l|. 
we have shown through VMC that the chiral d-wave RVB 
state on honeycomb lattice does not support topological 
degeneracy. However, since the calculation is done in the 
region A/x < v^, we would like to supplement it with 
the result in the region A/x > •s/Z- 

The overlap between states with different number of 
vison in both holes of the torus for the chiral d-wave 
RVB state is shown in Fig. 11. Here \^io) denotes the 
state with trapped vison in the hole surrounded by the x- 
circumference but no vison in the hole surrounded by the 
y-circumference. IvPoi)) denotes the state with trapped 
vison in the hole surrounded by the y-circumference but 
no vison in the hole surrounded by the x-circumference. 
The calculation is done at A/x = 5. The overlap is seen 
to oscillate with the lattice size. Such oscillation is caused 
by the complex nodal structure of the chiral d-wave RVB 
state, which move continuously in the momentum space 
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with A/x- However, it is clear that the overlap will ex- 
trapolate to finite value in the thermodynamic limit. 



-0.50- 




-0.54 -I , 1 , 1 , 1 , 1 , pi 

1 2 3 4 5 

FIG. 12: The variational energy per site for the Heisenberg 
model on honeycomb lattice calculated from the chiral d-wave 
RVB state. The calculation is done on a 16 x 16 x 2 lattice with 
periodic-antiperiodic boundary condition. The error bars are 
smaller than the size of the symbols. 



0.20 




0.05 - 



0.00- 1- l-l 

1 2 3 4 5 

FIG. 13: The spin structure factor at the antiferromagnetic 
ordering wave vector for the chiral d-wave RVB state on hon- 
eycomb lattice as a function A/x- The calculation is done 
on a 16 X 16 X 2 lattice with periodic-antiperiodic boundary 
condition. The error bars are smaller than the size of the sym- 
bols. The small non-analycity of the curve is checked to be 
finite size effect induced by the moving nodes in the momen- 
tum space. The red line denotes the mean field prediction, 
which is independent of A/x. 



E. Variational energy and spin structure factor 

The variational energy for the Heisenberg model on 
honeycomb lattice is shown in Fig. 12 as a function of 
A/x- As for the staggered flux phase on square lattice, 
the variational energy is optimized with a state interme- 
diate between uniform RVB state and the 7r-flux phase. 
Unlike the square lattice case, the difference in energy 
between the uniform RVB state and the optimized vari- 
ational state in this class of RVB state is seen to be 
small. As the exact ground state energy is estimated to 
be —0.55 J per site from exact diagonalizition studies[l8|. 
we see the chiral d-wave RVB state in fact represents a 
rather good variational state for the Heisenberg model 
on honeycomb lattice. 

The spin structure factor at the antiferromagnetic or- 
dering wave vector for the chiral d-wave RVB state is 
plotted in Fig. 13 as a function of A/x. As a compari- 
son we also plot the value calculated from the mean field 
ground state(= 4^), which is independent of the pairing 
strength. The spin structure factor is seen to be reduced 
at the uniform RVB state and the tt flux phase, both of 
which having a SU (2) gauge structure at the mean field 
level. This can be taken as a evidence that the SU{2) 
gauge fluctuation will reduce rather than enhance the 
spin correlation. 

We have also calculated the size dependence of the spin 
structure factor at A/x — 0, A/x — 1 and A/x = to 
see how the gauge structure affects the spin correlation 
in the thermodynamic limit. The results are shown in 
Fig. 14. We find in aU of the three states the spin struc- 
ture factor decay algebraically with the linear scale of 
the lattice, indicating no magnetic long range order in 
the thermodynamic limit. This is different from square 
lattice, where the uniform RVB state is found to posses 
antiferromagnetic long range order. A fit to the date with 
the formula S{Q) = a(l/i)^ shows that the exponent /3 
varies with the value of A/x • More specifically, the expo- 
nent at A/x = is found to be approximately 1.5. It is 
reduced to 1 at A/x = 1 and return to 1.3 at A/x = V^- 
Thus, a finite value for the flux will not only improve the 
local spin correlation, but also result in more long ranged 
spin correlation. 

V. CONCLUSION 

From the above discussion, we see the chiral d-wave 
RVB state on honeycomb lattice indeed stands as a nat- 
ural generalization of the staggered flux phase on square 
lattice. The two states shares the following properties. 

(1) Both states respect the full symmetry of the Heisen- 
berg model on respective lattices and have a positive def- 
inite wave function in the sense of the Marshall sign rule. 

(2) Both states evolves with A/x in the same manner. 
With the increase of A/x, both states evolve from the 
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A -3 -2 -1 



Ln(1/L) 

FIG. 14: The finite size scaling behavior of the spin structure 
factor of the chiral d-wave RVB state for A/x ~ (Dirac), 
A/x = 1 and A/x = \/2(7r-flux). 

uniform RVB state to the 7r-flux phase and then back 
to the uniform RVB state. The uniform RVB state and 
the TT-flux phase on both lattices have a promoted gauge 
symmetry of SU{2). 

(3) The introduction of the gauge flux on both lattice 
improves the short range spin correlation and both the 
staggered flux phase on square lattice and the chiral d- 
wave RVB state on honeycomb lattice are rather good 
variational description of the Heisenberg model on re- 
spective lattices. 

(4) As a result of the bipartite nature of the lattice and 
the resultant Marshall sign rule of the RVB state, both 
the staggered flux phase on square lattice and the chiral 
d-wave RVB state on honeycomb lattice do not support 
topological degeneracy, although the latter posses a Z2 
gauge structure at the saddle point level. 

However, there are also important differences between 
the staggered flux phase on square lattice and the chiral 
d-wave RVB state on honeycomb lattice. The following 
is a list of the main differences. 

(1) On square lattice, the intermediate state connect- 
ing the uniform RVB state and the 7r-flux phase has a 
U{1) gauge structure. However, on honeycomb lattice, 
the two states are connected by intermediate state with 
Z2 gauge structure. This difference originates from the 
frustrated nature of the dual lattice of honeycomb lat- 
tice. Related with this difference is the difference in the 
staggering pattern of the gauge flux on the two lattices. 
On square lattice, the gauge flux form a two sublattice 
collinear pattern. On the other hand, the gauge flux 
on the honeycomb lattice form a three sublattice non- 
collinear pattern. 

(2) As a result of the non-collinear nature of the gauge 
flux in the chiral d-wave state, the states with gauge flux 



$ below A/x = \/2 are not gauge equivalent to the states 
with gauge flux 27r — $ above A/x = a/2, except the spe- 
cial case of the uniform RVB state. Thus in principle the 
TT-flux phase on honeycomb lattice need not be a extreme 
of physical properties as functions of A/x, although we 
find both the variational energy and the spin structure 
factor do reach their extreme at the 7r-flux phase. On the 
other hand, on square lattice, the 7r-flux phase must be 
an extreme of physical properties as functions of A/x as 
protected by the $ 27r — $ symmetry. 

(3) On square lattice, the uniform RVB state posses 
antiferromagnetic long range order. However, on honey- 
comb lattice, no magnetic order is detected in the uni- 
form RVB state. This difference can be induced by the 
difference in their spinon dispersion. While the uniform 
RVB state on square lattice has a large and nested spinon 
Fermi surface, which is unstable toward antiferromag- 
netic ordering, the uniform RVB state on honeycomb lat- 
tice has an isotropic Dirac-type spinon dispersion, just 
as the TT-flux phase on square lattice. Indeed, the spin 
structure factor of the uniform RVB state on honeycomb 
lattice is found to show the same scaling behavior as the 
TT-flux phase on square lattice [2 2]. 

(4) Although the variational energy can be improved by 
introducing gauge flux on both square lattice and hon- 
eycomb lattice, the amount of improvement is quite dif- 
ferent. On square lattice, the introduction of staggered 
flux will lead to a 15% improvement in the variational 
energy of the Heisenberg model. On the other hand, the 
same energy gain on honeycomb lattice is less than 2%. 
This difference can also be induced by the difference in 
the spinon dispersion on the two lattices. 

As the staggered flux phase plays such a central role 
in our understanding of the high-Tc superconductivity in 
cuprates from the perspective of RVB theory, it is quite 
natural to expect that the chiral d-wave RVB state will 
play the similar role, if the doped honeycomb system in 
the strongly correlated regime does support some kind of 
superconductivity. After doping, the time reversal sym- 
metry will be broken and the state will become truly chi- 
ral. A full gap will also be opened in the doped system, 
as can be shown by including a chemical potential term 
in Eq.(13). There will be many interesting issues con- 
cerning this novel superconducting state, especially on 
its topological properties. However, as the condensation 
energy is much smaller than system on square lattice, 
the superconductivity on the honeycomb system should 
be weaker than high Tc cuprates. It is interesting to see 
if this state can be realized in experiment. 

This work is supported by NSFC Grant No. 10774187 
and National Basic Research Program of China No. 
2007CB925001 and No. 2010CB923004. The author is 
grateful for the discussions with Fan Yang and Cenke Xu. 
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APPENDIX 



In this appendix, we present the mean field spinon dis- 
persion of the chiral d-wave RVB state on honeycomb lat- 
tice for some typical values of A/x- The result is shown 
in Fig. 15. 
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FIG. 15: The evolution of the mean field spinon dispersion 
as a function of A/x- From top to bottom are the results for 
A/x = 0, A/x = 1, A/x = %/2, A/x = 2 and x/A = 0. 
The value of x and A are normalized so that x^ + A^ = 1. 
The lower branches(E]^) are plotted in the left panel and the 
upper branch (Ej^) are plotted in the right panel. 



